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ABSTRACT 

Recent observations suggest that long-duration 7-ray bursts (GRBs) and their 
afterglows are produced by highly relativistic jets emitted in core-collapse explosions. 
The pre-explosive ambient medium provides a natural test for the most likely pro- 
genitors of GRBs. Those stars that shed their envelopes most readily have short jet 
crossing times and are more likely to produce a GRB. We construct a simple computa- 
tional scheme to explore the expected contribution of the presupernova ejecta of single 
Wolf-Rayet (WR) stars to the circumstellar environment. Using detailed stellar tracks 
for the evolution of massive stars, we discuss the effects that the initial main sequence 
mass, metallicity, rotation and membership in a binary system have on the ambient 
medium. We extend the theory of GRB afterglows in winds to consider the effect of 
the relativistic fireball propagating through the WR ejecta. Specific predictions are 
made for the interaction of the relativistic blast wave with the density bumps that 
arise when the progenitor star rapidly loses a large fraction of its initial mass or when 
the ejected wind interacts with the external medium and decelerates. A re-brightening 
of the afterglow with a spectrum redder than the typical synchrotron spectrum (as 
seen in GRB 970508, GRB 980326 and GRB 000911) is predicted. We also calculate 
the luminosity of the reflected echo that arises when circumstellar material Compton- 
scatters the prompt radiation, and examine the spectral signatures expected from the 
interaction of the GRB afterglow with the ejected medium. 
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1 INTRODUCTION 

Recent evidence has given support to the idea that long- 
duration 7-ray bursts (GRBs) result from the cataclysmic 
collapse o f massive stars in very energetic supe rnova-like ex- 
plosions (MacFadyen, Woosley & Heger 2001) rather than 



collimation, the maximum isotropic 7-ray emission is raised 
to unprecedented values, implying a total energy of the or- 
der of 3 x 10 54 ergs (GRB 990123; Kulkarni et al. 1999). A 
beamed jet alleviates this energy requirement. The evidence 
is still preliminary, but collimation factors of Qj/in < 0.01 
have be en derived from the decline in some afterglo w light- 



the coa lescence of two compact objects, such as black holes curves (Kulkarni et al. 1999; Castro- Tirado et al. 1999 



or neut -on stars (Lattimer & Schramm 1976). The observed rison et al. 1999). Whether a jet is present or not, the energy 



Har 



distribution of optical afterglow luminosities with respect to 
their host galaxies suggests that some GRBs ma y be asso- 
ciated with star- forming regions ( Paczyhski 1998) ) . The spa- 
tial and temporal coinciden ce between the sup ernova (SN) 
1998bw and GRB 980425 (iGalama et al. 19981) hints that 



budget is compatible with current scenarios that involve the 
collapse of compact objects: the binding energy of the orbit- 
ing debris and the spin energy of the black h ole are the two 
main energy reservoirs available (Rees 1999). 



some GRBs may be linke d to ultra-bright t ype-Ibc SNe but 
this is still controversial (Pian et al. 2000). Some evidence 
for supe rnova-type emissio n has also been found in GRB 
980326 (|Bloom et al. 1999|), GRB 970228 ([Reichart 199E|) 



and GRB 000911 (Lazzati et al. 2001). If we believe that 
GRBs are associated with SN-like events, then the explo- 
sion is likely to be jet-ind uced if the baryon cont amination 
problem is to be avoided ( Rees fc Meszaros 1992 ). Without 



MacFadyen & Woosley (1999; hereafter MW99) have 
explored the evolution of rotating helium stars whose 
iron core collapse does not produce a successful tradi- 
tional neutrino-powered explosion. For massive stars (Mo i; 
35M0) with sufficient angular momentum, a rapidly accret- 
ing stellar mass black hole forms promptly at the centre of 
the collapsing star (Type I collapsars). Less rapidly accret- 
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ing black holes can also form over longer timescales due to 
the infall of stellar material that failed to escape during the 
initial SN (Type II collapsars). Prompt and delayed black 
hole formation can occur in stars with a range of radii de- 
pending on the evolutionary state of the massive progenitor, 
its metallicity and its membership in a binary system. If the 
star loses its hydrogen envelope along the way, and if the jet 
produced by the accretion maintains its energy for longer 
than it takes the jet to tunnel through the star, a GRB 
is likely to be produced (MW99). Otherwise, acceleration 
of the explosion debris to a sufficiently high Lorentz factor 
(> 10 2 , Meszaros & Rees 1997) is unlikely a nd an asym 



metric supernovae like SN 1998bw may result (MacFadyen 
Woosleffc Heger 200 1|) . 



In the collapsar model for GRBs, or in any other model 
involving a massive star, the key to obtaining relativistic 
motion is the escape of an energy-loaded fireball from the 
stellar environment. This is aided if the progenitor under- 
goes a Wolf-Rayet (WR) phase, which is characterized by a 
strong stellar wind that causes the star to lose enough of its 
outer layers for the surface hydrogen abundance to become 
minimal. The radius of a WR star is sufficiently small for 
the explosion energy to break out before the engine ceases to 
operate. Because of the intrinsic variations of mass-loss rates 
in WR evolution, the GRB blast wave expands into shells of 
varying gas density. The effects of the WR ejecta interacting 
with the interstellar medium (ISM) can be observed in the 
wind bubbles around some of these objects. The decelera- 
tion of a pre-SN wind by the pressure of the surrounding 
medium or the interaction of fast and slow winds (Luo & 
McCray 1991; Vikram & Balick 1998) could also create cir- 
cumstellar shells. 

It is an unavoidable consequence of a core collapse 
model for GRBs that the afterglow emission propagates 
in a dense stellar wind. Chevalier & Li (1999, 2000) de- 
scribed some features of the afterglow evolution in a steady, 
spherically symmetric wind with density oc r~ 2 . In this pa- 
per we study the effects of strong, non-steady WR stellar 
winds with large mass-loss variations and non-spherical ge- 
ometries. In particular we explore the impact of the pre- 
SN evolution of single high-mass star (Mo ;> 35M ) on 
the circumstellar environment. In § 2 we use detailed stellar 
evolution of WR stars to show how the progenitor initial 
mass (Mo) and metallicity Z affect the density profile of 
the ambient material. In § 3 we extend the theory of after- 
glows in winds to consider the propagation of a relativistic 
fireball through the dense environment expected at the end 
of the evolution of a WR star. We calculate the echo that 
arises when the progenitor circumstellar material Compton- 
scatters the prompt 7-ray radiation in § 4. Implications of 
recent observations, in particular for GRB 970508, GRB 
980326 and GRB 000911, and predictions for future ones 
are discussed in § 5. We examine the spectral line signatures 
expected from the interaction of a GRB afterglow with the 
chemically enriched environment in § 6. Our conclusions are 
presented in § 7. 



* The comparative modelling of single and binary systems will be 
presented in a subsequent paper (Dray et al. 2001 in preparation). 



2 EVOLUTIONARY MODELS OF MASSIVE 
STARS WITH HIGH MASS-LOSS RATES 

2.1 Model outline 

2.1.1 WR evolution 

WR stars have many extreme properties: they have lumi- 
nosities typically between 10 5 and 10 6 Lq and are very hot. 
This places them in the upper left corner of the HR diagram. 
Because of their dense stellar winds (10 _4 MQyr _1 ), the hy- 
drostatic surfaces of WR stars are mostly inaccessible to 
direct observations. This makes it difficult to associate stel- 
lar models with observed objects on the basis of their radii, 
surface temperatures or luminosities. The spectroscopic de- 
termination of the abundances of H, He, C, N and O (e.g. 
Smith & Hummer 1998) allows only certain classes of theo- 
retical stellar models. Massive stars are known to have stel- 
lar winds and mass-loss while on the main sequence and, 
although the amount of mass lost here may be small com- 
pared to that lost in later phases, it is relevant to the subse- 
quent evolution. The first evolutionary phase of a WR star 
is the nitrogen-line stage (WNL), which begins when CNO 
processed material is exposed at the stellar surface. This pro- 
gresses to the WNE stage during which no surface hydrogen 
is detectable. Once the hydrogen envelope is completely lost, 
the WR models strictly obey the mass-luminosity relation- 
ship for hydrogenless WR stars. This allows the details of 
their further evolution to be calculated by solving on ly the 
equation of global energy conservation (Langer 1987) for a 
given mass-loss rate, because the mass determines the lumi- 
nosity which in turn determines the rate of nuclear burning 
in the stellar core. After the helium envelope is shed, the WC 
stage occurs, in which strong carbon lines can be seen. Fi- 
nally, WO stars are formed, with high surface oxygen abun- 
dances. Subclassifications in the above categories are made 
according to line strength ratios (eg. Smith & Maeder 1991). 

In order to assess the simplest case, we shall confine 
the discussion in this section to non-rotating, non-magnetic 
single stars. We use detailed evolutionary tracks between 
the initial main sequence and the WO stage, with enhanced 
mass- loss rates based on the work of Meynet et al. (1994). 
We carry out the computational simulations with the stel- 
lar evolution code first developed by Eggleton (1971) and 
recently updated by Pols et al. (1995). Due to slight dif- 
ferences in the stellar evolution codes, i ncluding a differ- 
ent treatment of convective overshooting (Schroeder, Pols & 



Eggleton 1997), the stellar tracks computed by Meynet et 
al. (1994) are not identical to the ones calculated here (e.g. 
the minimum initial mass for WR star formation is slightly 
lower) but they are qualitatively the same. 

As described in Meynet et al. (1994), enhanced mass- 
loss rates were used during the phases when stellar winds 
are believed to be driven by radiation pressure, i.e. dur- 
ing the whole evolution with the exception of the WNE 
and WC/WO phases, which correspond to the stage of 
nearly pure helium envelopes and of He+C+O cores, re- 
spectively. Enormous mass-loss can occur at the entrance 
into these stages. Pre-WR mass-loss rates for the whole HR 
diagram were taken from the expression derived by de Jager 
et al. (1988) and scaled with metallicity as (M z /Mq) = 
(Z/Zq) 0,5 . For WNL stars a constant mass-loss rate was 
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adopted (Abbott & Conti 1987) and for later stages of WR 



evolution a mass-dependent mass-loss rate was assumed as 
specified by Langer (1989). In cases where the cumulative 
WNL mass-loss was enough to nearly evaporate the star we 
used the WNE rate instead. 



2.1.2 Wind-ISM interaction 

Winds from the most massive stars tend to have the greatest 
effect on the ISM because the mass-loss rates are large and 
the winds are very fast and carry large momentum fluxes. 
Moreover, massive stars tend to lie in associations, so there 
is the cumulative influence of the winds of many stars. In 
describing the effects of stellar winds on the surrounding 
medium we used the basic theory for the wind-ISM interac- 
tion develop by Castor, McCray & Weaver (1975). We as- 
sume that the ISM has a simple structure, with no clumps, 
no magnetic fields, no multiple ionization phases and no 
nearby stars. During the evolution of a wind-driven circum- 
stellar shell the system has a four-zone structure (analogous 
to that of a supernova shell; Woltjer 1972; Cox 1972). From 
the inside to the outside these zones are: (I) a supersonic 
stellar wind with density p(r) = M/47rr 2 Woo; (II) a hot, al- 
most isobaric region consisting of shocked stellar wind mixed 
with a small fraction of the swept-up interstellar gas; (III) a 
thin, dense, cold shell containing most of the swept-up inter- 
stellar gas; (IV) ambient interstellar gas of number density 
no. 

The wind initially expands unopposed into the ISM 
with a velocity of about Woo, the escape velocity at the 
sonic point (for the purpose of this analysis we assume 
Woo = yj2GM* / R*, because the sonic point is reached after 
a few stellar radii). The free expansion phase is considered 
to be terminated at a time idee, when the swept-up mass 
of the interstellar medium is comparable to the mass in the 
wind (Woltjer 1972; Castor et al. 1975; McCray 1983). The 
mass lost by the star is Mt& ec and the swept-up mass is 
-j-(i>oo tdec) 3 nom P Hc, where p c ~ 2 in a helium gas. These 

two masses are equal when tdcc = \J 3Af/(47ro§onom p /i c ), 
which is about 100 years for a typical WR wind expand- 
ing into an homogeneous ISM^. During this time the wind 
bubble has reached a radius of 



o o w mW i>l/2 —1/2 — 1/2 

r dcc = 8.8 x 10 M_' s n 0fi < v^ ^ cm, 



(1) 



where M is the mass loss rate in units of solar masses per 
year, no is the density of the surrounding medium in units 
of cm -3 , Woo the wind velocity in units of km s _1 and we 
adopt the convention Q — 10 x Q x . 

In a dense (no ~ 10 5 cm~ 3 ) molecular cloud a wind- 
driven bubble will radiate its energy and stall in a short 
time, causing the ejected wind to slow rapidly at a much 
smaller radius (Shull 1982). The same may be true if a neigh- 
bouring star is present - for example, a very close binary in 



which the WR wind is strongly slowed down by the OB 
star radiation pressure before it collides with the OB wind. 
-Rob < tob < fter, where roB is the distance from the OB 
star (with radius -Rob) to the region of the stellar wind colli- 
sion and r tC r ~ 5 Rob is the distance from the OB star cen- 
tre at which the OB wind reaches terminal velocity^]. Many 
WR stars with absorption lines appear to be single; how- 
ever, the fraction of visible close WR + OB binaries seems 



to be around 35 per cent (Chiosi & Maeder 1986) 



When the free expansion phase (I) has ended, the wind 
encounters an inward facing shock. Kinetic energy is de- 
posited in the shocked wind region in the form of heat, 



3 m p p c 2 
Tshock = -T7 — ; — (Aw) 
lb k 



An 



L4xlQJ( 1 nnv -i ) K ' 
100 km s 

where Aw is the relative speed of the material approach- 
ing the shock. This creates a temperature of about 10 7 K in 
the shocked wind region. During phase II, the material is 
so hot that it causes the contact surface to move outward 
more slowly than it would in a freely expanding wind. The 
ISM that enters the outward facing shock is heated to a tem- 
perature below 10 K, emission of line radiation becomes the 
dominant cooling process and the swept-up gas cools quickly 
to temperatures of about 10 4 K that can be maintained by 
the radiation field of the star. The duration of the adiabatic 
expansion phase can thus be estimated by finding the time 
it takes the expanding gas to cool from T s h oc k ~ 10 7 K to 
10 6 K. Using equation (2), we find that a change in temper- 
ature from 10 7 K to 10 6 K corresponds to a change in jump 
velocity by a factor of \/T0. This change in jump velocity 
represents a phase 11:1 age ratio of about 6 (McCray 1983). 
Thus, the age of the adiabatic phase is less than about 1000 
years. 

The mass of the swept-up material is much larger than 
that in the hot wind and, because it is cool, it lies in a com- 
pressed region. Phase III persists for as long as the star is 
able to sustain a powerful wind. The dominant energy loss 
of region II is work against the compressed region III (see 
Castor et al. 1975). The compressed region III expands be- 
cause its gas pressure is higher than that of the surrounding 
ISM. Therefore, the expansion is described by the momen- 
tum equation, 

d 
dt 



[M s (t)v(t)]='i7vr 2 (t)P i , 



(3) 



where Pi is the internal pressure of the compressed region, 
assuming that most of the swept-up interstellar mass re- 
mains in the thin shell. Ms{t) is the mass of the shell of 
swept-up material, given by Ms(t) = (4/3)7rr 3 (t)po. Pi is 
determined by the gas pressure of the high temperature gas 
in region II. The wind material that enters the backward 
facing shock is hot, but the material that enters the for- 
ward facing shock is cool. The cooled swept-up material is 
driven outward by the high gas pressure of the hot bub- 
ble. The stellar wind adds energy to region II at a rate 



t The free expansion phase takes place at the early stages of 
the evolution of the hot star and occupies a minimal fraction of 
its lifetime. During this time both M and Vca are approximately 
constant (see Figs, la and lc). 



t If a very close binary is present, the high density region will 
probably be reduced to a small covering factor surrounding the 
collapsing star. Note that a beamed fireball would introduce an 
additional geometrical factor. 
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Figure 1. Stellar winds from different evolutionary WR models. The state of the ambient medium of a WR star at the end of its life 
(d) is determined primarily by the time evolution (until core collapse) of the mass-loss rate (a) and the escape velocity of the wind (c). 
Panel (b) shows the mass of the WR star as a function of time. Stellar models were computed for initial masses of 40 {solid line), 50 
{dotted line) and 60 {dash-dotted line) Mq. The external medium particle density of the wind ejected by the progenitor is best modelled 
(at r < 

r dca sce equation 1) by a power-law n — At ^. A spherical shell with high column density arises when the ejected material 
accumulates at the radius at which the free expansion phase of the ejected wind terminates. We have set the density of the ambient 
material before the WR phase to 1cm -3 . 



iw(t) = ^A/(t)n^ (t). The internal energy in the bubble is 
given by the product of the energy per mass of the ma- 
terial, (3/2)nkT/ pi — (3/2)Pi/pi, and the total mass of 
the bubble, (4/3)-7rr pi. Since the total internal energy of 
the bubble comes from the energy of the wind, we find 
Pi = L w {t)/[2nr 3 {t)] (Castor et al. 1975; McCray 1983). 
The expansion of the bubble during the adiabatic phase can 
be found numerically by using this expression in the mo- 
mentum equation^]. The bubble could continue to expand 
until stalled by the pressure of the ISM. However, before 
that happens, the star will have produced a supernova type 
explosion which re-pressurizes the bubble. 



2.1.3 Wind-wind interaction 

During the different mass-loss stages within a WR lifetime, 
the wind power of the star increases or decreases (as shown 
in Fig. la) with time. When the escape velocity of the wind 



S If the wind power L w is roughly constant for a period of time, 
t, one can write Pi = L w t/(27rr 3 ). The resulting solution of equa- 
tion (3) gives r{t) oc i 3 / 5 . This shows that the shell expands more 
slowly than would a freely expanding wind. 



increases, the fast wind collides with the early ejecta before 
reaching r^cc- The resultant shell will be pushed outward by 
the central star wind and retarded by the early ejected wind, 
quickly reaching a constant velocity, but increasing in mass. 
The basic theory for the wind-wind interaction was devel- 
oped by Kwok, Purton & Fitzgerald (1978). The resulting 
expansion law for the shell follows from the momentum bal- 
ance between the two components (denoted as i and j). If 
Mij is the mass of the resulting shell when it is at a radial 
distance rij{t), then 



r« (t) M 3j 

Mij{t) = I -^-dr + 



Mi 



dr. 



(4) 



where Vj < Vi. If Vij(t) is the velocity of the shell, then, 
assuming a completely inelastic collision, the equation of 
motion may be written 



. d Vi 



Mi 



M, 



M iS (t)^- = ^-[v i -v ij ] 4 -^-[v i 



(5) 



Numerical integration of equation (5), with a substitution 
for Mij{t) from equation (4), gives the resulting expansion 
law for the shell. The thickness of the shell Ar^- may be 
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steepness of the hydrogen profile above the H-burning shell, 
the WNL stage is rather short in low initial mass stars; its 
duration increases towards higher masses and it may last for 
t he whole po st-main-sequence phase for the highest masses 
(Langer 1987). Due to the existence of a H-burning shell and 
also because hydrogen has a high radiative opacity and a low 
mean molecular weight, the radii of WNL stars are consider- 
ably larger (and their surface temperatures correspondingly 
smaller) than WNE and WC stars. The escape velocity of 
the wind therefore decreases during this phase, as shown in 
Fig. lc. During the early phase of WR evolution, the pro- 
genitor stars lose a large fraction of their initial mass quite 
rapidly, as shown in Fig. lb. We have studied the interaction 
of the strong stellar winds with the surrounding interstellar 
gas and have found that a dense circumstellar shell devel- 
ops at the radius at which the free expansion phase of the 
ejected wind terminates, as illustrated in Fig. Id. For a given 
metallicity, this overdense region is more extensive for higher 
initial masses. 



I n most early discus s ions ( Meszaros, Rees fc Wijers 



Figure L2 Ttio cffcrfg r.f HifWnnt initial maaaoa nn thn gtriif- 



lOOit IChcvalicr fc Li 1999j; Madau, Blandford fc Rees 200C 
Panaitescu fc Kumar 2000D , the progenitor star was ex- 



ture of the ambient medium. The density profile can be fit by a 
power law n(r) = Ar~ 2 for r < r'dcc- The value of the constant 
A as a function of Mo is shown in panel (a) for the stellar models 
computed in § 2.1. A was scaled to A=3.0xl0 35 A* cm -1 . The 
radial location of the highest density shell as a function of Mo is 
shown in panel (b) for different values of no. The free expansion of 
the ejected wind is terminated when the swept-up mass becomes 
comparable to the mass of the wind. The ejected mass then accu- 
mulates at such a radius (r(j cc ). This overdense region lies closer 
to the progenitor for WR stars originating from low initial mass 
stars. All the evolved WR stars have lost their hydrogen envelope 
and are left with a bare helium core (see panel c). 



found by requiring its internal pressure to balance the pres- 
sure from the wind (see Kwok et al. 1978). 



2.2 The effects of mass-loss during the evolution 
of WR stars 

An increase of the mass-loss rate has profound consequences 
for the modelling of the WR population. At a given metallic- 
ity, the minimum initial mass for the formation of a WR star 
(lower for higher M), the duration of the WR stage (greater 
for higher M) , the times spent in the different WR subtypes 
and the surface composition during these phases are all very 
sensitive to the mass-loss rates. A detailed discussion is given 
in Maeder (1991). Here we emphasise some interesting con- 
tributions to the structure of the ambient medium that arise 
when the ejected winds interact with the ISM during the 
WR lifetime (about 10 5 yr). 

Using the physical ingredients mentioned in § 2.1, we 
have computed the evolution of stars with initial main- 
sequence masses, Mo, between 10 and 150 Mq and of solar 
metallicity, ^0=0.02. Fig. la shows the time evolution (until 
core collapse) of the mass-loss rate for Mo = 40, 50, and 60 
Mq. The progenitor star loses a large fraction of its initial 
mass quite rapidly on entering the WNL phase (the high- 
est constant mass- loss rate interval in Fig. la). The main 
reason is the non-negligible contribution of the H-burning 
shell to the total luminosity. Depending somewhat on the 



pected to be surrounded by a substantial 
n(r) — Ar~ a 



(6) 



medium at the end of its life (with a = 2 for a wind ejected at 
a constant speed). For the stellar models computed here, we 
found no significant deviations from an r~ 2 density gradient 
as a result of the small mass-loss variations from the pro- 
genitor star prior to core collapse. We calculated power-law 
indices, s, with values ranging from 1.9 to 2.1. By running 
stellar models with a variety of initial masses, we found that 
the constant A in the density profile is not a strong function 
of Mo, as shown in Fig. 2a. The constant A in equation (6) 
was scaled to A = 3.0 x 10 35 A* cm 1 as in Chevalier & Li 
(2000). The highest deviation from the extrapolated density 
profile is mainly provoked by the interaction between the 
winds and the interstellar medium. As time progresses the 
amount of material sweptup increases and the momentum 
from the wind cannot drive the shell to such high speeds. 
The ejected mass then accumulates at a corresponding ra- 
dius, creating an overdense region. We found the radial lo- 
cation of the overdense region or bump to be smaller for 
WR stars that originate from low initial mass stars (see Fig. 
2b). The radius at which the free expansion ends is a func- 
tion of the product M t ^v^^ri ~ , so that one may fear 
that large variations of no can be masked or mimicked by 
much smaller variations of Mo (or M 1 ^ 2 v^^ 2 ; see Fig. 2b). 
It should be remarked that the plausible range of no (0.1 - 
10 5 cm -3 ) is large compared to the plausible range of Mo 
(Mo > 35 Mq if GRBs originate from rapidly rotating mas- 
sive stars with large helium cores, Fig. 2c). 



2.3 The role of metallicity 

Metallicity Z influences the stellar evolution of massive stars 
mainly through bound-free and line opacities. The main 
source of opacity in the interior is electron scattering, which 
is independent of metallicity. Therefore, metallicity has little 
direct structural effect in massive stars. Bound-free and line 
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Figure 3. The effect of metallicity on stellar winds. An 85 Mq main sequence star was evolved with three different metallicities: Z = 0.5 
Zq (solid line), Z = 0.05 Zq (dotted line), Z = 0.015 Zq (dash-dotted line). Low metallicity decreases the mass-loss rate (a), increases 
the mass of the helium core (b) and keeps the radius of the star small, reducing the escape velocity of the stellar wind (c). Also, low 
metallicity favours high mass- loss rates at the end of the WR phase, increasing the density of the ambient medium (d) at r < r^ cc . A 
spherical shell with high column density arises when the ejected wind interacts with the ISM (no = 1 cm -3 ) and decelerates. In the case 
of the 85 Mq main sequence star with Z = 0.015 Zq a thin shell is present at r < r^ cc due to the substantial change in mass-loss prior 
to core collapse. 



opacities, on the other hand, are important in the outer lay- 
ers of massive stars. In this way, metallicity influences th e 
mass-loss rates by stellar winds (Chiosi & Maeder 1986). 



Low metallicity keeps the radius of the star smaller and 
also reduces the mass-loss. Both effects inhibit the loss of 
angular momentum from the evolving star (MW99). For a 
given mass-loss rate, the lower the metallicity, the higher the 
stellar mass for WR formation. Furthermore, low metallicity 
raises the threshold for the removal of the hydrogen envelope 
by stellar winds, thus increasing the mass of the heaviest he- 
lium core and favouring black hole formation (MW99). 

We have computed the evolution of stellar models for 
an initial main sequence mass of 85 Mq at three different 
metallicities: Z = 0.5 Zq, Z = 0.05 Z G and Z = 0.015 
Zq (see Fig. 3). As discussed above, low metallicity reduces 



the mass- loss for most of the WR lifetime. However, low 



metallicity favours high mass-loss rates at the end of the 
WR evolution (see Fig. 3a). This increases the density of 
the ambient medium at small radii prior to core collapse 
(see Fig. 3d). Nevertheless, the overall amount of ejected 
material that accumulates at r^cc increases with metallicity, 
as shown in Fig. 3d. 

By running stellar models with a variety of main- 



sequence initial masses and metallicities, we found that the 
normalization constant of the density profile, A,, is not a 
strong function of Z, as shown in Fig. 4a. The highest de- 
viation from the extrapolated profile occurs at progressively 
smaller radii for WR stars originating in low metallicity en- 
vironments (see Fig. 4b). In such environments, the forma- 
tion of massive helium stars, desired for GRB formation, is 
favoured (see Fig. 4c). 

The numbers of WR stars in galaxies of different metal- 
licities are an important test of the stellar models at various 
Z and of the values of the final stellar masses. The theoret- 
ical predictions of the ratios of WR/O-stars, WC/WR and 
WC/WN are in very good agr eement with observations of 
the Local Group (Maeder 1991). This supports t he adopted 
dependence of the mass-loss rates on metallicity ( Meynet et 



al. 1994). At high Z (in inner galactic locations), gas opaci- 



ties are larger in the outer stellar layers and so more momen- 
tum is transferred by radiation pressure, mass-loss is more 
intense and thus GRB formation is disfavoured. The strong 
dependence of massive star evolution on metallicity suggests 
that GRB progenitors in inner galactic locations may be in- 
trinsically different from those in outer ones. Furthermore, 
this dependence on Z implies that GRB characteristics could 
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Figure 4. Stellar winds from WR models with different mctallic- 
ity. All the evolved WR stars have lost their hydrogen envelope. 
The value of the normalization constant A„ as a function of Mo is 
shown in panel (a) for different values of Z. The highest density 
shell lies closer to the progenitor for WR stars originating from 
lower Z (see panel b). Low metallicity raises the threshold for 
removal of the hydrogen envelope by stellar winds, increasing the 
mass of the helium core (panel c). 



be strongly affected by redshift (MW99). 



2.4 The effects of WR rotation 

If long GRBs are SN-like events, the explosion is likely to 
be the outcome of black hole formation through the collapse 
of a rapidly rotating massive star (WM99; Heger, Langer & 
Woosley 2000). Rotation may affect WR evolution in sev- 
eral ways (Maeder & Meynet 2000, and references therein). 
First, in the case of fast rotation, the star may enter the 
WR phase while still burning hydrogen in its core and thus 
spend more time there. Second, rotation (through its ef- 
fects on both the mass-loss rates and mixing) may favour 
the formation of WR stars with lower initial mass. Typi- 
cally the minimum mass for WR formation is 35-40 Mq for 
non-rotating model s. It decreases to about 25 Mq for initial 
Kot = 300 km s" 1 ( [Maeder fc Meynet 1987| ). Finally, higher 
rotation-velocities lead to longer WR lifetimes and higher 
mass-loss rates. Thus, rotation could alleviate the need to 
enhance the mass-loss rates in order to reproduce the ob- 
served WR/O-star number ratio as originally proposed by 
Meynet et al. (1994). 



Direct attempts to measure the rotation velocity of WR 
stars have been performed only for a few cases: Massey 
(1980) obtained wsini ~ 500 km s" 1 for WR138 (see also 



Koenigsberger 1990); Massey & Conti (1981) measured 
vsini ~ 150 — 200 km s" 1 for WR3. Also, indirect evidence 
has been found for the existence of some axisymmetric fea- 
tures around WR stars. For example, Harries et al. (1998) 
suggest that about 15 per cent of WR stars have anisotropic 
winds caused by equatorial density enhancements produced 
by rapid rotation. High rotation velocities for iron cores 
favour the collapsar model for GRBs (see MW99). The con- 
tribution of the WR mass-loss to the surrounding medium 
along the GRB beam, which is likely to be close to the ro- 
tation axis of the star, may be smaller than calculated here 
for the spherical wind case. 



3 WIND INTERACTION MODELS FOR GRB 
AFTERGLOWS 

One can understand the dynamics of the afterglows of GRBs 
in a fairly simple manner, independent of any uncertainties 
about the progenitor systems, from the relativistic gener- 
alization of the method used to describe supernova rem- 
nants. The basic model for GRB afterglow hydrodynamics 
is a relati vistic blast wave expand ing into the surrounding 
medium (VIeszaros & Rees 1997). The interaction of the 



outer shell with the external medium is described by the adi- 
abatic Blandford-McKee (1976, hereafter BM) self-similar 
solution. The scaling laws that are appropriate for the burst 
interaction with a medium with particle density n oc r~ s 
have been described by Meszaros et al. (1998), Chevalier & 
Li (1999,2000), and Panaitescu & Kumar (2000). We aim 
to examine the specific predictions for the interaction of the 
relativistic blast wave with the density bumps that arise 
when ejected stellar material decelerates against the exter- 
nal medium (see § 2). This collision would be an important 
contribution to the observed afterglow if it were not to take 
place too late in the fireball evolution. We outline the rele- 
vant physical parameters that determine when this collision 
occurs in § 3.1. We discuss briefly the situation when the 
burst is collimated in § 3.2. We analyse the shock conditions 
when the shell collides with an enhanced density medium 
in § 3.3, and estimate the synchrotron emission that results 
and its observational consequences in § 3.4. 



3.1 Hydrodynamics of a relativistic shell 

For an adiabatic ultra-relativistic blast wave, the (isotropic 
equivalent) total energy is 



E 



87rArV 



2„3- 



17 -4s 



(7) 



where F is the bulk Lorentz factor of the shock front and 
r is the observed radius near the line of sight (BM). We 
assume the burst to be collimated with an initial half-angle 
8 larger than 20°, and that lateral expansion is negligible 
during the relativistic phase. A distant observer receives a 
photon emitted along the line of sig ht towards the fireba ll 
centre at a time t = r/4(4 - s)r 2 c ( |Ctievalier fc Li 2000|) , 
and so 



{4- s) (17 -4s) Et 



2-kAc 



(8) 
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Before the collision with the high density shell, the shock 
front is expected to propagate through an n(r) — Ar~ s 
wind. In a spherically symmetric wind ejected at a constant 
speed the density drops as s=2. Using the stellar models in 
§ 2, we found 0.1 < A, < 10 for WR stars surrounded by 
an s m 2 medium at the end of their life (A t — A / 3.0 x 10 
35 cm -1 ). For these stellar parameters, a re-brightening of 
the afterglow as a consequence of the collision of the shock 
front with the high density spherical shell is observed at a 
time 



tda 



0.25 



E 



1 Q52 er gg 



U0 17 c 



At 



(9) 



after the burst. Here tday is the observer time measured in 
days. If the collision takes place at a sufficiently late phase 
of the afterglow the system is adiabatic. This is probably 
valid from about half an hour after the burst (Granot, Piran 
& Sari 1998). The outer shell catches up with the thin high 
density shell, which is at rest at a radius r s h n. The shell col- 
lision produces two new shock waves: a forward shock that 
moves into the thin shell and a reverse shock that propa- 
gates into the relativistic ejecta. The calculation of shocks 
in the collision of two cold shells resembles the calculation of 



energy emitted from internal shocks ( Kumar & Piran 2000 ) 



3.2 Collimation of 7-ray bursts 

The energy flux from a GRB may be collimated, as wit- 
nessed in the jets from active galactic nuclei and some bi- 
nary black hole sources. Some models explicitly include this 
collimation (MW99). As long as the jet is highly relativis- 
tic, the observed features are reproduced by spherical mod- 
els, but as the shocked jet slows down there are clear con- 
sequences for the observed afterglow (Rhoads 1997): the 



edge of the jet becomes visible when V = (2(9) ~ and the 
jet is able to expand laterally, leading to an exponential 
slowing of the forward shock front. In the case of a colli- 
mated GRB going off in a medium with density decreasing 
as r -2 , jet effects are expected to become important at a 



time tjet « 2(0/O.l) i E 52 AZ 1 days (|Chevalier fe Li 2000D 



Kumar & Panaitescu (2000), have recently pointed out, how- 
ever, that a fireball propagating into an s — 2 medium 
should show little evidence for light-curve steepening due 
to edge and sideways expansion. This could explain the lack 
of breaks in the afterglow lightcurves of GRB 980326 and 
GRB 980519. In a collimated outflow the sharpest break in 
the lightcurve occurs in a uniform low density medium and is 
associated with the edge of the jet approaching the relativis- 
tic beaming cone. In wind scenarios, where no clear break 
in the lightcurve is predicted, jets can perhaps b e detected 
by time depende nt measurements of polarization ( Kumar & 
Panaitescu 2000j) . 



For a collimated outflow expanding in an r 3 medium, 
the observed radius near the line of sight is proportional to 
{Et/e 2 Ac) 1/(A - s) . Thus, for a fixed value of E, the impact of 
the relativistic shell with a density bump would take place 
earlier in the evolution of the fireball if the burst is colli- 
mated, causing a re-brightening of the afterglow to appear 
at earlier times in the observer's frame. 



3.3 Shock conditions 

Consider a relativistic shell with a Lorentz factor r\ moving 
into a cold external medium (EM). The interaction is de- 
scribed by two shocks: a forward shock propagating into the 
EM and a reverse shock propagating into the shell. There 
are four regions separated by the two shocks: the EM (1), 
the shocked EM (2), the shocked shell material (3) and the 
unshocked shell material (4). The EM is at rest relative to 
the observer. From the jump conditions for the shocks and 
the equalities of pressure and velocity along the contact dis- 
continuity we can estimate the Lorentz factor 7, the number 
density n and the energy density e in the different shocked 
regions as functions of the three quantities ni, 714 and r\ 
(BM). The unshocked EM has 71 = 1, while the unshocked 
shell material moves at the original coasting velocity of the 
particles, 74 = 77. The afterglow emission is determined by 
the deceleration time of the shell. The emitting region moves 
toward the observer with a Lorentz factor P = 72 =73, as- 
suming that the shocked material (either region 2 or 3) emits 
radiation on a timescale shorter than the hydrodynamical 
time. 

Ther e are two limits in wh ich there is a simple analytical 
solution ( |5ari fe Piran 1995 ). If the shell density is high, 
W4 3> T] 2 ni, then the reverse shock is Newtonian (or mildly 
relativistic). In this case the energy conversion takes place in 
the forward shock, and the collision is too weak to slow down 
the shell efficiently, V ~ rj. If the density of the expanding 
shell is low, 714 <C rfrii, then the reverse shock is relativistic, 
the shell material decelerates considerably and 



(10) 



The energy density in the shocked regions satisfies e = e2 = 
eg — 4r 2 nim p c 2 . Similar relations hold for the reverse shock 



( [gari fe Piran 1999 ). In any realistic situation the EM is 
probably inhomogeneous, as in the stellar models in § 2. 

Consider a density jump by a factor a over a distance 
I. The forward shock propagates into the EM with a density 
ni as before, and when it reaches the position where the EM 
density is an\ a new shock wave is reflected. This shock is 
reflected again off the shell. Sari & Piran (1995) showed that 
the reflection time is « 1/iAca 1 ^ 2 ). After these reflections the 
Lorentz factor and hydrodynamical properties of the system 
are as if the EM were homogeneous with a density am . The 
corresponding observed timescale due to this inhomogeneity 
of the EM is t w l/(a 1/2 F) (for l£r). 



3.4 Synchrotron emission 

The synchrotron spectrum from relativistic electrons that 
are continuously accelerated into a power-law energy dis- 
tribution comprises four power-law segments, separated by 
three critical frequencies, the self-absorption frequency (f S a), 
the cooling freq uency (u c ) and the characteri s tic synchrotron 
frequency (z/ m ) (Sari, Piran & Narayan 1998; Meszaros, Ree: 



fe Wijers 1998 ). The spectrum and light-curve of an after- 
glow are determined by the time evolution of these frequen- 
cies, which in turn depend on the hydrodynamical evolution 
of the fireball. 

The break frequencies can be calculated if the energy 
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Figure 5. The effect of the impact of a relativistic shell with the 
density discontinuity on the optical afterglow. The shock front ex- 
pands within the s = 2 stellar wind until it reaches the high den- 
sity shell at a distance r^eii ~ 4x 10 17 cm. The radial profile of the 
high-density shell is modelled as the inner bump shown in Fig. Id 
for the 40 Mq case. Collisional models were computed for a high 
density enhancement a=100 (solid line) and a low density en- 
hancement a=10 (dotted line). At the time of the collision the rel- 
ativistic shell Lorentz factor is V ~ 3 for £52=0. 5, A*=l, £ e =0.1, 
es=1.0, and p=3. The afterglow emission is calculated in the adi- 
abatic regime. The collision model takes into account the fireball 
geometrical curvature when calculating the photon arrival time 
and relativistic boosting. Given the radius of collision r s j, e n S> I, 
the observed variability timescale is At Q b s ~ f s hell/( c r 2 ) > 20 
days. Note that T <C 3 after the collision. 

distribution of the injected electrons and the strength of the 
magnetic field are both known. The distribution of the in- 
jected electrons is assumed to be a power law of index — p, 
above a minimum Lorentz factor ji m e c (m p / m e )Y . The en- 
ergy carried by the electrons is a fraction (p — — 2) of 
the total internal energy e e . The turbulently generated mag- 
netic field, assumed to be amplified up to a fraction eb by the 
processes in the shocked region and not determined directly 
by the field in the WR star, is B « y/ 32ircSg 2 nip n J^T 
( |5ari, Piran fe Narayan 1996 ). For an adiabatic blast wave, 



the corresponding observer peak frequency is 



B 



2-p 



(ii) 



The ratio of v m , a the observed peak frequency from the rela- 
tivistic shell after the forward shock has traversed a density 
jump of a factor a over a distance I, and the frequency in 
the absence of such a jump v m , is 



rl(am) 1 / 2 



-1/2 



1/2 



(12) 



Therefore, we expect this emission to be a very significant 
contribution to the short-wavelength flux. The spectrum at 



low frequencies (y <C u m ) scales as u 1 ^ 3 ; at high frequen- 
cies (v 2> Vm) it scales as i/ - ^ -1 -*/ 2 . At late times, when 
the collision has run its course, the shells have merged back 

to a BM solution and the observed flux is proportional to 
£(p+i)/[2(4-s)] £ or n ^ r -s ^ j n aDsence f a d en _ 

sity enhancement, the flux would have been proportional 
to E^ +1)/[2{A - s)] (Kumar & Piran 200d). Therefore, the in- 



crease in the observed emission from the forward shock due 
to the density jump is approximately (E a / E) (p+1)/l2(4 ~ a)] « 
Q ,(p+ 1 )/l 4 ( 4 - s )] ( neglecting the enhanced energy losses from 
the shocks that may arise during the collision with the high- 
density thin shell. We should expect, however, a smooth 
transition from one solution to the other at most wave- 
lengths. This is not the case at frequencies for which there 
is significant emission from the reverse shock, from which 
the total emission is smaller, and at a lower frequency by 
typically a factor of F. At these frequencies the impact of 
this additional emission is likely to be significant. 

We conclude that the impact of the relativistic shell 
with the density discontinuity (see Figs. Id and 3d) pro- 
duces a bump in the lightcurve that can be very dramatic 
at low frequencies. To illustrate this, we consider the scal- 
ing laws described by Chevalier & Li (1999) for an adiabatic 
blast wave in an s = 2 medium. Before the collision, the 
synchrotron emission frequency of the lowest energy elec- 



trons IS Vr, 



5 x 10 1 



-i t s, 



1/2 ,-,1/2.-3/2 



-l rj 52 1 



day 



Hz and the 



.1/2 £,1/2.-1/2 
''day 



mjy 



flux at this frequency is F Vm = 20 A, \ ~\E 52 *d 
(These expressions assume z = 1 in a flat Universe with 
Ho — 65kms _1 Mpc -1 ). For electrons with a power-law dis- 
tribution p, the flux above v m is F v = F„ m (v /v m )~^ p ~ 1 ^ 2 oc 
£-(3p-i)/4 Q ver ^ ne typical time of observation of a GRB 
afterglow, the shock front expands within the stellar wind 
until it reaches the high-density bump at a distance of about 
4 x 10 17 cm. The impact of the relativistic shell with the 
density enhancement produces a clear feature in the after- 
glow lightcurve. Fig. 5 shows the effect for two different val- 
ues of the bump density. The lightcurve of the enhanced 
emission arising from the interaction of the forward shock 
with the density bump is calculated with the contribution of 
the curvature to the temporal structure described in Feni- 
more, Madras & Nayakshin (1996). Given a typical radius of 
collision r s heil ^ fcheil, the observed variability timescale is 
Ai bs ~ 7" s hcii/(F 2 c). The maximum enhancement of the ob- 
served emission with respect to the extrapolated afterglow 
lightcurve is approximately a 1 / 2 and the characteristic syn- 
chrotron frequency is lower than that in the absence of the 
collision. 



4 COMPTON ECHOES 

If massive stars are the progenitors of GRBs, then the hard 
photon pulse will propagate in a pre-burst stellar wind. The 
GRB is assumed to have a broken power-law spectrum, 
eE t oc e (e < 250 keV) and e -0 25 (e > 250 keV). Cir- 
cumstellar material Compton sca tters the prompt radiation 



and g ives rise to a reflection echo ( Madau, Blandford fc Ree, 1 



200C| ). If E = J E^ndedQ, is the total energy emitted by the 
burst, where E e n is the energy emitted per unit energy e 
and unit solid angle along the 9 direction, then the equiv- 
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alent isotropic luminosity of the Compton echo inferred by 
a distant observer is 



(13) 



where r is the distance from the site of the burst, da/dD. is 
the differential Klein-Nishina cross section for unpolarised 
incident radiation, ot is the Thomson cross section, e is 
the energy of the scattered photon and n e the local electron 
density. The equal-arrival time scattering material lies on 
the paraboloid determined by 



ct 



1 — cos 



(14) 



where 9 is the angle between the line of sight and the direc- 
tion of the reflecting gas as seen by the burst. 

Madau et al. (2000) derived analytical expressions for 
the echo luminosity lightcurve in the Thomson regime for 
a range of possible emission geometries and ambient gas 
distributions. In particular, they discussed the effect of an 
infinitesimal high-density spherical shell in the surround- 
ing medium. For isotropic emission, the echo luminosity de- 
creases as the equal-arrival time paraboloid sweeps up the 
shell, reaching a minimum when t = r s h e ii/c, and then in- 
creases until it reaches the back of the shell. For a collimated 
burst, the echo lightcurve is the sum of two delta functions, 
as the result of the interaction of the approaching and reced- 
ing jet with the spherical shell (see equation 10 in Madau et 
al. 2000). The two spikes of emission are seen separated by 
an interval 2r s h e ii cos 8/c, where 6 is the angle between the 
line of sight and the approaching beam (i.e., the main burst 
emission is not seen by the observer) . 

Fig. 6 shows the reflected echo of a two-sided GRB jet 
propagating through the dense environment expected at the 
end of the evolution of a 40 Mg Wolf-Rayet star with Zq 
(see Fig. Id). We calculated the luminosity of the reflected 
echo by integrating equation (13) with the full Klein-Nishina 
cross section. The GRB is assumed to radiate a total of 10 53 
erg, with each jet having equal strength. On the equal-arrival 
time paraboloid, the receding beam is reflected by denser gas 
that is closer to the source, and so its contribution domi- 
nates the echo at all energies where the scatter occurs in the 
Thomson regime. This is clearly illustrated by the second 
spike of the echo caused by the interaction of the reced- 
ing jet with the density bump in the surrounding medium. 
The spectral energy distribution mirrors the prompt burst 
at low energies, but is much steeper above a few hundred 
keV (Madau et al. 2000). The hard X-ray flash from a pulse 
propagating into the environment expected at the end of the 
evolution of a 40 Mq star could be detectable by Swift out 
to z ~ 0.1. It is instructive to look at the contribution of the 
density bump to the luminosity of the reflected echo. The 
time interval between the the two spikes, 2r a heii cos 8/c, can 
be used to determine the angle between the line of sight and 
the approaching beam. 



RECENT OBSERVATIONS AND FUTURE 
PREDICTIONS 




10 
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Figure 6. The Compton echo of a GRB. The primary burst is as- 
sumed to be a two-sided collimated pulse of energy 10 52 ergs s _1 
and duration 10 s, propagating at an angle 8 to the line of sight, 
and invisible to the observer. The GRB jet propagates through 
the dense environment expected at the end of the evolution of a 
40 M Wolf-Rayet star with Zq (see Fig. Id). 



5.1 Afterglow Sources 

The variety of observed af terglows, while compatible wit h 
relativistic fireball models (Meszaros, Rees & Wijers 1998), 



poses challenges for interpretation. Some of the key ques- 
tions include the effect of the external medium, the possible 
anisotropy of afterglow emission and the radiative efficiency. 
In this paper, we have considered the effects of the surround- 
ing medium due to massive WR stellar progenitors on GRBs. 
We found that a high-density bump arises when the ejected 
wind interacts with the external medium and decelerates, or 
when the progenitor star rapidly loses a large fraction of its 
initial mass. The impact between the relativistic shell and 
this density bump could be observed as a re-brightening^ 
of the afterglow with a spectrum typically redder than the 
synchrotron afterglow spectrum, as seen in GRB 980326 and 
GRB 970228. These observations have been explained as an 
underlying supernova outshining the afterglow some days af- 
ter the burst event. However, this interpretation should be 
regarded as tentative, in consideration of the fact that only 
R band measurements and a low signal-to-noise spectrum 
were obtained for GRB 980326. In the case of GRB 970228 
no spectroscopic measurement was made at the time of the 
suggested re-brightening, and the lightcurve was unevenly 
sampled in different filters. 

> Note that a re-brightening can also be produced by shock re- 
freshment due to delayed energy injection by an extremely long- 
lived central engine. In this case, however, the spectrum would be 
bluer than a typical GRB afterglow (see Panaitescu, Meszaros & 
Rees 1998 for the case of GRB 970508). 
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Figure 7. Modelling of the multi-filter lightcurves of the after- 
glow of GRB 000911. From top to bottom the normalised photo- 
metric measurements are plotted, scaled by factors of 10 2 , 10 1 , 
1, 10 — 1 , 10~ 2 for clarity. The dashed lines show the best fit for 
a collisional model obtained by computing the interaction of a 
rclativistic blast wave with the WR ejecta. 



red spectrum F v oc v~ 5 (compared with F v oc v~ 1A in the 
extrapolated afterglow) can be better reproduced if there 
is a substantial contribution to the emission from the re- 
verse shock at lower frequencies. We have fitted this model 
to the simultaneous multi-filter data (see Fig. 7), obtain- 
ing a reduced \ 2 intermediate between the fit with a simple 
external shock synchrotron component plus a galaxy model 
(x 2 ~ 1-3) and the fit with the additional SN component 
(X 2 ~ 1.1, see Lazzati et al. 2001). 

5.2 Internal shock models 

In the internal shocks scenario for GRBs, the actual 7-ray 
temporal profile is the outcome of the complex dynamics 
of the ejecta. A growing consensus is that a central site re- 
leases energy in the form of a wind or multiple shells over a 
period of time commensurate with the observed duration of 
GRBs (Rees & Meszaros 1994). Each subpeak in the 7-ray 
lightcurve is the result of a separate explosive event at the 
central site. In this early phase, the time-scale of the burst 
and its overall structure follow, to a large extent, the tem- 
poral behaviour of the source. In contrast, the subsequent 
afterglow emerges from the shocked regions of the external 
medium where the relativistic flow is slowed down and the 
inner engine cannot be seen directly. A concern was raised 
that internal shocks without deceleration were rather ineffi- 
cient (Kobayashi, Piran & Sari 1997) converting at most 25 
per cent of the bulk motion energy into radiation. Since the 
afterglows can only account for a few percent of the radiated 
energy, it is unclear where most of the energy goes. 



5.1.1 GRB 970228 & GRB 980326 

In order to explain the re-brightening of GRB 980326, we 
would require a large density jump with a size comparable 
to the variability timescale, about 10r lightdays w 1 17 cm, 
and located at a distance 20r 2 c days ~ 10 18 cm (Bloom 



et al. 1999). These conditions are clearly compatible with 
the density jumps calculated in § 2. A very red spectrum 
F v oc v~ 2 S (compared with F v oc v~ 10 in the extrapolated 
afterglow) can be produced if there is some contribution to 
the emission from the reverse shock at lower frequencies. In 
the case of GRB 970228 at z = 0.695, the isotropic lumi- 
nosity is 6 x 10 51 erg s _1 . The temporal bump detected 10 
days after the explosion can be explained by the impact of 
a relativistic shell expanding within an n(r) = A„r~ 2 stel- 
lar wind with the density discontinuity located at a distance 
5 x 10 17 cm (see equation 9). These values are in agreement 
with the overdense regions calculated in § 2. 



5.1.2 GRB 000911 

Using detailed stellar tracks for the evolution of massive 
stars, as described in § 2, we found that the re-brightening 
of GRB 0009111 at z = 1.1 (Lazzati et al. 2001), can be 
explained by the interaction of a relativistic blast wave 
(e e =0.1, £b=0.1, and p=2.2) expanding into the ambient 
medium expected at the end of the life of a 40 Mq (7 Mq 
core) WR star evolving with Z=Zq. The shock front ex- 
pands within a n(r) = 1.5 x 10 35 cm -1 r -2 stellar wind until 
it reaches the density enhancement (see Fig. Id). The very 



In a dense molecular cloud the density of the surround- 
ing medium, before the free expansion phase, can be as large 
as n « 10 4 - 10 5 cm" 3 (Shull 1982; Sanders, Scoville & 
Solomon 1985) This is not unlikely because massive stars 
tend to cluster and these stars do not live long enough to 
travel far from their place of birth. In this environment the 
free expansion phase of the ejected wind, which is termi- 
nated when the swept-up mass becomes comparable to the 
mass in the wind, ends at a radius of 10 15 — 10 16 cm (see 
equation 1). The ejected mass then accumulates at such a 
radius, creating an overdense region. Substantial internal 
energy can be converted into radiation at this radius when 
the shells responsible for the 7-ray emission run into this re- 
gion. This is equivalent to the model proposed by Fenimore 
& Ramirez- Ruiz (2001), with the dense region acting as the 
decelerating shell. No substantial deceleration is required to 
convert up to 80 per cent of the bulk-motion energy dur- 
ing the GRB. Indeed, a recent analysis of 387 pulses in 28 
BATSE GRBs shows that the most intense pulses within a 
burst have nearly identical widths throughout the burst, but 
that the weak pulses tend to become wid er as the burst pro- 
gresses ( Ramirez- Ruiz & Fenimore 2000 ) . This effect is seen 



in internal shock simulations when a high density thin shell 
at a radius of 10 15 — 10 16 cm is included. Furthermore, the 
detection of an opt ical flash associated with GRB 990123 
(Akerlof et al. 1999), with a peak magnitude of 9 in the V 



band 50 s after the initial burst (successfully described by 
the reverse shock model by Sari & Piran 1999), can also be 
explained as a result of the impact of two shells t hat collide 
with a rela tive Lorentz factor of the order of 2 (Kumar & 
Piran 200C). Such collisions are expected in internal-shock 
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models if an ejected shell catches up with a decelerating 
shell or when a shell runs into a dense region at a radius of 



which Fc absorption edges and Kq features would be present 



(Mcszaros fc Rees 1998 A.mati et al. 2000) 



6 THE 7- RAY BURST ENVIRONMENT 

The task of finding useful progenitor diagnostics is simpli- 
fied if the pre-burst evolution leads to a significantly en- 
hanced gas density in the immediate neighbourhood of the 
burst. The detection of spectral signatures associated with 
the GRB environment would provide important clues about 
the triggering mechanism and the progenitor. Stars inter- 
act with the surrounding interstellar medium, both through 
their ionising radiation and through mass, momentum and 
energy transfer in their winds. Mass-loss leads to recycling 
of matter into the interstellar medium, often with chemical 
enrichment. Mass- loss is a particularly significant effect in 
the evolution of massive stars; in particular, for stars that 
enter WR stages. WR stars follow all or part of the sequence 
WNL, WNE, WC and WO, which corresponds to a progres- 
sion in the exposure of nuclear products: CNO equilibrium 
with H present; CNO equilibrium without H; early visibil- 
ity of the products of the 3a reaction; and then a grow- 
ing (C+0)/He ratio respectively. Mass-loss drastically in- 
fluences stellar yields. At low Z there is a high production 
of a-nuclei, while at higher Z most of the He and C pro- 
duced is ejected in stellar winds and escapes further nuclear 
processing. 

The environment could also be metal-enriched, al- 
though solar abundances are probably more likely. The 
amount of metals in the environment can be larger than 
the stellar yields if a supernova explosion occurs before the 
GRB, perh aps producing significant Fe Kq and K-edge lu- 
minosities ( Lazzati, Campana fc Ghisellini 1999| ) . Because of 
the very high luminosities and hard initial 7-ray spectrum, 
all the Fe is fully ionised, and as it cools, the strong Fe Kq 
and K-edge appear initially in absorption, and la ter as re- 
combination emission features (Weth et al. 2000). Guided 
by the report of Fe line detectio ns that peak after 1 day in 



GRB 97050 8 and GRB 970828 (Piro et al. 1999 



Yoshida 



et al. 1999), we are led to consider a high density shell. A 



physical requirement is that the distance ctV(t) , reached 
after 1 day by the afterglow shock, is less than the shell ra- 
dius. Notice that in order to reproduce the observed Fe Kq 
equivalent widths, a shell located at few 10 16 cm is required 
with a mass of Fe ~ 2.5 x 1O _4 M0, in a total shell mass 
of 1 M e (Weth et al. 2000). Such enhanced density shells 
can be present in the environment of a progenitor star if the 
density of the surrounding medium, before the WR phase, 
is of order 10 5 cm -3 . The ejected mass accumulates at this 
radius, creating a high-density region 10 9 -10 11 cm -3 , for a 
50 Mq main-sequence star. Without substantial decelera- 
tion of the wind ejecta, such high density shells are difficult 
to form. Emission lines are hard to detect, unless the event 
occurs in an exceptionally dense or metal-rich environment 
(Piro et al. 2000). A massive star environment could be dis- 
tinguished by the presence of a significant flux of Fe Ka and 
probably H Ly-a emission lines, reprocessed by a moderately 
Thompson thick envelope. For a more massive thick enve- 
lope, a significant reflected component would be expected, in 



7 CONCLUSIONS 

Recent observations strongly suggest that at least some 
GRBs are related to massive supernovae. Quantitative in- 
sight into the formation of GRBs is hindered by the lack 
of detailed core-collapse calculations. The ground-breaking 
work of MW99 suggests that GRBs are more likely to oc- 
cur in stars that have lost their hydrogen envelope. Stars 
with less radiative mass-loss retain a hydrogen envelope in 
which a poorly collimated jet that loses its energy before 
breaking through the stellar surface is likely to arise. Highly 
relativistic jets will not escape red supergiants with radii 
> 10 13 cm. By contrast, a focused low-entropy jet that has 
broken free of its stellar cocoon is likely to arise from a WR 
progenitor. Our detailed stellar tracks for the evolution of a 
WR star show the diverse effects that initial main-sequence 
mass, metallicity, and rotation can have on the surrounding 
medium. Pre-explosion mass loss provides a natural medium 
to generate GRB afterglows, and allows specific predictions 
about their temporal structure to be made. The presence of 
density bumps in the nearby ambient environment are in- 
herent to the evolution of WR stars. The impact between 
the forward shock and these high-density regions should be 
observed as a re-brightening of the afterglow with a typically 
redder spectrum. Interestingly, the density ring in the wind 
profile lies closer to the progenitor for WR stars with low 
initial mass and low metallicity. This characteristic offers a 
direct observational test of which stars are likely to produce 
a GRB. 

The total energy observed in 7-rays from GRBs whose 
redshift has been determined is diverse. One appealing as- 
pect of a massive star progenitor is that the great variety of 
stellar parameters can probably explain this diversity. Given 
the need for a large helium core mass in progenitors, the 
burst formation may be favoured not only by rapid rotation 
but also by low metallicity. Larger mass helium cores might 
have more energetic jets, but it is unclear whether they can 
be expected to be accelerated to large Lorentz factors. Many 
massive stars may produce supernovae by forming neutron 
stars in spherically symmetric explosions, but some may fail 
neutrino energy deposition, forming a black hole in the cen- 
tre of the star and possibly a GRB. One expects various 
outcomes ranging from GRBs with large energies and dura- 
tions, to asymmetric, energetic supernovae with weak GRBs. 
The medium surrounding a GRB would provide a natural 
test to distinguish between different stellar explosions. 
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